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Rigidity and stability of cold dark solid universe model. 
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The large scale dynamics of the universe appears to be dominated by a "dark energy" constituent 
with negative pressure to density ratio w = P/p, which could be stable if sufficiently rigid, but not 
if purely fluid. It was suggested by Bucher and Spergel that such a cosmological solid might be 
constituted by a cold (static) distribution of cosmic strings with w — —1/3, or membranes with 
the observationally more favoured value w = —2/3, but it was not shown that the rigidity actually 
would be sufficient for stability. For cases in which the defect lattice is formed from even junctions, 
it is found that the rigidity to density ratio will be given by p/p = 4/15 in both the string and 
membrane cases, and it is confirmed that this is indeed sufficient for stabilisation with respect to 
sufficiently small perturbations. 

PACS numbers: 95.35.+d, 98.80.-k 
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INTRODUCTION 

It was pointed out by Bucher and Spergel Q that 
an underlying dark energy component consisting of a 
frozen network of topological defects, having an approxi- 
mately uniform (unclustered) density distribution p and 
a strongly negative value of its pressure P, could ac- 
count for many of the observed features of the universe 
on a large scale. In particular, the cosmic acceleration 
suggested by measurements of type la supernovae. It 
was suggested that, due a symmetry breaking phase 
transition at a cosmological temperature of a few hun- 
dred KeV or less, such a constituent might be provided 
by a suitable distribution of membranes (domain walls), 
whose averaged tension would be able to provide a nega- 
tive pressure, provided the positive pressure contribution 
from kinetic effects is not too large, i.e. provided the dis- 
tribution is effectively cold enough to be treated as an 
approximately static "frozen" state. 

For ordinary Nambu-Goto type cosmic strings, Bucher 
and Spergel noted Q, Q that in a static isotropic distri- 
bution the pressure to density ratio, w = P/p (in units 
with unit light speed) would be given by w = —1/3. 
Moreover, for ordinary Dirac type membranes (such as 
simple domain walls) in a cellular lattice the value would 
be w = —2/3, which is fairly close to what is consid- 
ered 0,0 to be observationally favoured. 

The optimism of earlier 0] and some more re- 
cent 0, 0, 01 discussions of such a proposal has been 
challenged from various points of view. One form of ob- 
jection concerns conceivable observable consequences 
derived from assumptions that are themselves open to 
question. However, it would seem that the most funda- 
mental problem is the absence, so far, of any convincing 
suggestion as to how an effectively "hot", wiggly, disor- 
dered structure - such as would naturally emerge from 



a Kibble-type spontaneous symmetry breaking - could 
be damped in such a way as to "freeze" as an effectively 
"cold" lattice of straight string or flat wall segments, as 
is required for the pressure to become effectively neg- 
ative. The positive pressure contribution from wiggles 
and waves will of course tend to be removed naturally by 
standard damping mechanisms, but the problem is that 
such processes will usually do so by removing the 
segments involved (so that the mean pressure to density 
ratio remains positive) , and will not necessarily tend to- 
ward a regular lattice. 

The present work does not attempt to address this ba- 
sic problem of finding a plausible mechanism for setting 
up a cold string or membrane system of the postulated 
kind in the first place. Our purpose here is just to provide 
a better understanding of the subsequent mechanical be- 
haviour of such a system - assuming it has already been 
somehow created - and in particular to investigate con- 
ditions under which such a negative pressure system can 
be stable. 

For a purely fluid medium, negative pressure would 
of course entail local instability on a short timescale, as 
symptomised by an imaginary value for the sound speed 
c s given by 



dP 
dp~ 



(1) 



In response to this implied objection, Bucher and Spergel 
pointed out that a medium constituted of cosmic strings 
or membranes would behave like an elastic solid, rather 
than a perfect fluid, which can be stable even when the 
pressure is negative, provided the rigidity is sufficient. 

The aim of the present work is to justify this claim, 
at least for membrane systems which contain only even 
junctions. However our calculation does not apply to 
wall systems with odd junctions, which we expect to be 




FIG. 1: Illustration of two dimensional lattice of the odd type, 
to which the present analysis does not apply. Each of the 
three colours corresponds to a different vacuum state. This 
hexagonal lattice is provided by a complex scalar field model 
of the form © with M = 3. 

much less rigid, so that they are likely to be, typically, 
too unstable for relevance as a dominant dark energy 
constituent. (If the walls were sufficiently massive, due 
to formation by symmetry breaking at very high energy, 
systems with odd junctions might nevertheless have been 
relevant to the primeval structure formation mechanism 
discussed in refs. [UE3)- O ur quantitative estimates will 
be based on the supposition that it is a sufficiently good 
approximation to treat the system as effectively isotropic 
on a macroscopic scale (meaning large compared with 
the mesoscopic lattice separation scale) but this simpli- 
fying condition does not seem essential for our qualita- 
tive conclusions, which we expect to remain valid even 
after allowance for small deviations from macroscopic 
anisotropy. 

DISTINCTION BETWEEN "EVEN" AND "ODD" 
TYPE SYSTEMS 

The calculation we present here will specifically be con- 
cerned with membrane systems that can be described as 
even-type, meaning systems where the number of walls 
intersecting at each (string-like) junction is even. Alter- 
natively, if the number of walls meeting at each junction 
is odd, this will described as of odd-type. Mixed systems 
containing both odd- and even-type junctions might arise 
in sufficiently complicated field models, but systems of 
that kind have not been put forward to date and will not 
be considered here 

Simple odd- and even-type systems can be obtained in 
two dimensions from complex scalar field models with a 
potential 

V{4>} = A|$ M -ry M | 2 , (2) 
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FIG. 2: Analogue of Fig.^with four differently coloured vac- 
uum states obtained from (|5J by setting M = 4. This case 
provides a lattice of the even type to which the present anal- 
ysis applies. 

where M, A and rj are constants. There are m minima of 
this potential at $ = 77exp[im/A/] for m = 0, 1, .., M — 
1. For M = 3 one can construct a hexagonal lattice 
as illustrated in Fig. ^ which is clearly of the odd-type, 
whereas for M = 4 an even lattice as illustrated in Fig. [21 
is possible. We note that the stability of such systems, 
which can be derived from a super-potential [lf|, has 
been investigated in ref. [l2^ . 

Our analysis will be based on the postulate that the 
affine displacements completely describe the modifica- 
tion of the equilibrium distribution at a mesoscopic level 
(meaning large compared with the membrane thickness, 
but small compared with the lattice separation scale), 
that is, we will assume that the local equilibrium of the 
system will preserved by arbitrary affine space trans- 
formations. We describe this as the affine equilibrium 
preservation postulate. 

This postulate is applicable to lattices of the even- type, 
for which the flat membrane sheets cross straight through 
each other (as in an X-form) at an angle that can freely 
adjusted without affecting the condition of equilibrium. 
However, it can be expected to fail for odd-type systems: 
for example at a (Y-form) junction between 3 equivalent 
membrane segments in equilibrium the intersection angle 
must always be 27r/3, a condition that will not be pre- 
served by a deformation that is affine on a mesoscopic 
scale. In such a case, a macroscopically affine defor- 
mation would therefor entail non affine (rigidity lower- 
ing) adjustments at the mesoscopic level. The effect of 
a locally affine (shear wave) transformation on the two 
dimensional complex scalar field model with M = 4 is 
illustrated in Fig. |31 

More complicated, and possibly stable, systems of both 
odd and even type type can be obtained from perturbed 
O(N) models with potential energy density given by N 
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FIG. 3: The effect of a (verticaffy propagating) shear wave 
of simple mesoscopically affine type on the even type lattice 
illustrated in Fig|5] It is clear that the local equilibrium of 
each of the junctions is maintained under this locally affine 
deformation. This would not be possible for the hexagonal 
lattice illustrated in Fig. for which a shear deformation 
could be affine only at a macroscopically averaged scale. 



real scalar fields $i , and mass scales rj , £ , in the form 



(3) 



where A and £ are dimensionless parameters such that 
£ > —A , NA > —£. The exactly O(N) symmetric model 
is obtained by taking £ = , A > 0, while the simplest 
example of an even-type system is provided by the de- 
coupled limit, A = , £ > 0, for which the domain walls 
simply pass through each other without interaction. 

Less trivial possibilities include the systems with £ = 
considered by Kubotani and collaborators 0, 0], who 
specifically envisaged odd-type systems - with triply in- 
tersecting (Y-form) junctions - obtained for N=3 by tak- 
ing £ < 0. In this case there are 6 vacua which can 
be thought of as being at the centre of the faces of 
a cube. The alternative, of relevance for the present 
work, is the even case - with quadruple (X-form) junc- 
tions [f| between walls separating vacuum domains where 
$2 = ( Ar; 2 + £f)/(N\ + £) for i = l,..,N - which can 
be obtained (for any value of N) by taking the symmetry 
breaking parameter to be positive, £ > . For N=3 there 
are 8 vacua at the vertices of a cube. 

Field models like this engender lattices of periodic or 
quasi-periodic type, with walls grouped in mutually par- 
allel families, the prototype example being that of the 
cubic lattice for N=3, in which there are just 3 families 
arranged so as to be mutually orthogonal. The question 
of whether disordered (glass-type) configurations can be 
obtained from other plausible models remains unclear. 
Although its averaged stress tensor will be isotropic, the 
cubic (N=3) case (see Fig. 0J) will be characterised by 
an elasticity tensor that distinguishes between directions 




FIG. 4: Illustration of a 2 dimensional section through a 
frozen cubic lattice of "even" type, having (differently shaded) 
domains bounded by flat walls with X form junctions, ob- 
tained for system with eightfold vacuum provided by broken 
O(N) model with e > and N=3. Such a configuration has an 
averaged stress tensor that will already be exactly isotropic, 
but to get a very highly isotropic elasticity tensor a rather 
larger value of N would be needed. 



parallel to the planes of the walls and directions that can 
deviate from these planes by as much as ir/4. For larger 
values of N, and correspondingly larger numbers of fami- 
lies of mutually parallel walls, the possible deviations will 
get progressively smaller so that, like the stress, the en- 
suing elasticity properties will become highly isotropic, 
as postulated in the analysis below. 



LOCAL STABILITY CONDITION 

For a system of the even-type just described, a lattice 
of flat membrane states continuing straight through each 
other - as if without interaction - at simple crossover 
junctions will evidently provide a local equilibrium con- 
figuration in a flat space background. Such an equilib- 
rium will clearly be stable at the mesoscopic scale char- 
acterised by the lattice cells, because the unbounded - 
effectively non-interacting - membrane sheets are each 
individually stable with respect to local perturbations, 
provided their amplitude is small compared with the lat- 
tice spacing, so as to prevent distinct vacuum domains of 
the same type from colliding. 

Assuming that - at a macroscopic scale large com- 
pared with the cell spacing - such a lattice interacts with 
other (relatively lightweight) constituents of the cosmic 
background sufficiently to form an effectively coherent 
medium (like a gas of point particles on scales large com- 
pared with the mean free path) one would expect that 
this medium would also be stable, which means having 
enough effective rigidity to provide real values for the 
relevant squared velocities of propagation. 
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To verify that, for an approximately isotropic unper- 
turbed configuration, this will indeed be so, we need the 
formula for the effect of rigidity on the relevant propaga- 
tion velocity. This was first studied [l3j when relativistic 
elasticity theory was originally dev elop ed for application 
to the solid crust of a neutron star [l4j , but the same for- 
malism applies equally well here. As in the terrestrially 
familiar non-relativistic limit, there will be transversely 
polarised - shake type - modes with propagation speed 
c ± , as well as a longitudinally polarised - sound type - 
mode with propagation speed c,, . 

It was shown |l3| that the latter would exceed the or- 
dinary fluid sound speed value Q by an amount propor- 
tional to the rigidity modulus p according to the formula 



2 , 4 2 



in which the squared shake mode speed is given by 



c 2 = 



/' 



(4) 



(5) 



- P + P 

For stabilisation of the medium, that is to ensure that 
c 2 will be positive even if c 2 is not, the necessary and 
sufficient condition is therefore 



(6) 



where is the bulk modulus, as defined in the relativistic 
case ^1 0| by the formula 



dP 

P = {p + P) — 
dp 



(7) 



This reduces to the familiar form = p dP/dp in the non- 
relativistic limit \P\ <C p. The requirement JSJ simplifies 
to p > — 37P/4 in the case of a polytropic equation of 
state, defined in terms of a conserved number density n 
and constants 7, K, and to by p = en 1 + mn which 
implies P = (7 — I) kti 1 . 

The models considered here correspond to the zero- 
mass limit, to = 0, in this polytropic class, and the bulk 
modulus will be given by 



= wjp, 



w 



7-1 



(8) 



Since we have 7 = 2/3 in the string case, and 7 = 1/3 
in the membrane case, it follows that the bulk modulus 
will have a negative, and thus destabilising, value that 
will be given by the same formula, 



(9) 



in both string and membrane cases. It can be seen that 
this result will still be valid for a non-polytropic mixture 
with p = k% ra 1 / 3 + ft 2 n 2 / 3 . Therefore, either separately 
and also for a mixed system of strings and membranes, 



the stability criterion JBJ will reduce to a requirement of 
the same simple form, 
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EVALUATION OF RIGIDITY 



(10) 



The new result provided by the present work for the 
even-type lattices described above is the explicit evalua- 
tion of the rigidity p, and the confirmation that the ratio 
p/p - which turns out to be the same for strings as for 
membranes - actually will satisfy the stability condition 

en. 

The required rigidity, or shear, modulus p can be de- 
fined, for an initially isotropic state, by expressing the 
static - quadratic order - change in the mass-energy den- 
sity p due to an infinitesimal volume conserving space 
displacement as 0,0,0] 



Ap = p e i:j e l - 



(11) 



where e,j is the infinitesimal strain tensor 



For an 

affine displacement given in Cartesian coordinates by 
x % 1— > x l = x l + e*jX? , it is defined as the symmetric 
part of the deformation matrix: = (e^ + £jj)/2. 

In order to avoid an extra energy variation contribu- 
tion due to the pressure, we impose the restriction that 
the displacement must conserve volume. This requires 
that the Jacobean determinant \dxi/dxj | should be unity 
which implies that the trace e\, like the variation Ap, 
will vanish to linear though not higher - order in e. 
An obvious way to do this is to take the symmetric (zero 
curl) pure shear deformation given in terms of coordi- 
nates x 1 = x, x 2 = y, x 3 — z by x == (I + e)x, y = y, 
z = z/(l + e), which when substituted in lfTT|) gives 



Ap = 2/i e 2 



(12) 



to second order in e. The value of p can be read out from 
this when the variation Ap has been worked out. 

However, the required value of p can be obtained in 
a computationally more convenient way by considering 
the effect on p of the asymmetric (non-zero curl) simple 
shear deformation given by 



x = x + 2 e z , 



y = y, 



(13) 



which is of the kind produced dynamically by a transverse 
mode with propagation speed in the z direction, given by 
|[SJ. It can be checked that a simple shear deformation 
such as this will provide a second order expression of the 
same form l|f 2|) to the one obtained for the pure shear 
deformation considered above. 

To evaluate the density variation Ap produced - as 
the left hand side of (|f 2f) - by the action of the simple 
volume conserving transformation l|f 3|) on a distribution 
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of randomly oriented strings, let us consider the case of 
a sample string segment from the origin to a point with 
coordinates given by the components, {x,y, 1} say of a 
vector £, which will be mapped to a vector £ + A£ with 
coordinates {x + 2 e, y, 1}. The length I = (x 2 + y 2 + 1) 1 / 2 
of the segment will undergo a fractional change given, to 
quadratic order in e, by 



A£ 
T 



2xe 2{£ 2 -x 2 )e 2 



P < £ 4 < 14 > 

Since the mass-energy of a static Nambu Goto string seg- 
ment is obtained simply by multiplying its tension, T 
say, the formula (|14|1 immediately provides the fractional 
change in the contribution to the mass density p from 
strings with the chosen direction, as given in polar coor- 
dinates by x = tan 9 cos <f>, y = tan 9 sin <p. 

Since the model is based on the postulate that the 
strings have a random isotropic distribution, their net 
effect can be computed by just taking the spherical aver- 
age of ()14J) over the polar coordinate values 9, 4>, which 
gives a formula in which the first order part cancels out, 
leaving just the quadratic term, 



Ap 
P 



At 
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(15) 



To evaluate the effect of the same simple shear defor- 
mation 1(13|) on an isotropic distribution of plane mem- 
branes, consider a parallelogram shaped sample seg- 
ment spanned by two vectors, £ 1 , with components 
{XjjljO}, which will not change, and £ 2 , with com- 
ponents {x 2 ,0, 1}, which will be deformed to {(x 2 + 
2e), 0, 1}. The membrane analogue of the string formula 
(|15fl for the effect on the density of an isotropic distribu- 
tion is 



Ap _ /AA\ 



(16) 



where A is the sample area. This is given by the 
magnitude of the vector product A = £ 1 x £ 2 , with 
components {l,x 1 ,—x 2 } in the direction normal to the 
membrane. The deformation 1(13(1 will map this area- 
calibrated normal vector to A + A A, with components 
{1, x 1 , — (x 2 +2 e)}. Since this transformation differs from 
that of the string segment considered above only by a per- 
mutation of the roles of the x and z axes, the final result, 
after spherical averaging, will be the same. 

For the macroscopically isotropic and mesoscopically 
even-type systems under consideration, this reasoning 
provides a theorem to the effect that - like the bulk mod- 
ulus (3 as given by © - the rigidity modulus p of the 
isotropic cosmic membrane distribution will be identical 
to that of an isotropic cosmic string distribution with the 
same density p. Thus, by comparing 1(15(1 with 1)12(1 . it 
can be seen that the value of the rigidity is 
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P- 



(17) 



not only in the string case, but also in the membrane case 
and for a mixed system of strings and membranes. 



CONCLUSIONS 

It is evident that the value given by 1)17(1 does indeed 
satisfy the stability condition 1(10(1 . This confirms the 
intrinsic coherence of the solid universe model subject to 
the availability of an underlying field theory of a suitable 
kind, meaning one providing a membrane system of even- 
type, for example, J3J with £ > 0. However, it remains 
unclear how such a system could have come into existence 
in the first place. 

In a solid of this kind, as a consequence of and J5J, 
we obtain real values for the speeds of longitudinal as 
well as transverse modes. For the latter, (JjJJ with 117|) 
gives c 2 = 4/157, so the (transverse) shake wave speed 



will be given by 
7 = 2/3 => c 



2/5, 7=1/3 



4/5. 



for the string and membrane cases respectively The cor- 
responding formula obtained from Q for the longitudinal 
wave speed will take the form 



w + 



16 

457 ' 



(18) 



which will indeed be consistent with the stability crite- 
rion (J5J), since it clearly provides a real- valued propaga- 
tion speed for any positive value of the polytropic index, 

7 = w + 1 . 

In particular it can be seen that the longitudinal propa- 
gation speed will be given by 

w=-l/3 => c 2 = l/5, w = -2/3 =>• c 2 = 2/5, 

in the string and membrane cases respectively, the latter 
being the one of greatest cosmological interest. 

The present analysis is based on the cold limit in which 
thermal 0] or other kinetic excitations of the world 
sheets are neglected. If significant, such effects would 
tend to reduce the effective tension, giving slower 
perturbation speeds and also making it harder to explain 
the apparent [3j cosmic acceleration. Such effects would 
also tend to undermine the stability of the system, which 
is valid (as shown above) only for perturbation ampli- 
tudes that are small compared with the relevant lattice 
separation scale. The fundamental difficulty with this 
kind of scenario is the problem of finding a natural mech- 
anism for getting rid of such excitation effects so as to 
actually obtain a frozen lattice of the required kind. 

Assuming that the membrane domination scenario can 
be achieved, and that it could be obtained by some spon- 
taneous symmetry breaking process, a first attempt was 
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made Q at evaluating the energy scale, r\ say, charac- 
terising the tension, T m 77 , of the required membranes 
and the cosmological temperature, 6 « 77, of their pre- 
sumed formation. It relied on the assumption that the 
relevant correlation lengthscale at the epoch of formation 
would be the maximum permitted by causality. In a re- 
cent reconsideration of this issue by one of the present 
authors it has be shown that, although it would be 
more realistic to suppose the correlation length to be 
very much smaller than the causal limit, nevertheless the 
effect of the appropriate correction on the estimate for 
77 is relatively moderate: the likely range is merely re- 
duced from the order of a hundred KeV to a few KeV. 
This does not change the qualitatively important con- 
sequence that the transition should have occurred at a 
relatively recent cosmological era, after pair annihilation 
and Big-Bang Nucleosynthesis, though before the epoch 
of recombination when the universe became transparent. 
The corollary is that the present lattice separation scale 
of the membrane lattice will be comparable 0, 01- at 
most, with interstellar rather than intergalactic separa- 
tion lengthscales. Thus, on the presumption that the rel- 
evant fields interact very weakly with ordinary particles 
such as electrons, and particularly photons, so that the 
walls will not be directly visible, Bucher and Spcrgcl's 
provisional conclusion that their scenario poses no hier- 
archy problem remains valid, in the sense that the range 
of conceivable values for 77 overlaps with the mass range 
for already known elementary particles. 

Despite this attractive feature, the plausibility of such 
a scenario is undermined by the absence, as remarked 
above, of any mechanism for getting from a randomly 
excited initial state to a regular frozen membrane lattice 
of the required kind, as exemplified by the case depicted 
in Figure 1. In view of the positive results of our study, 
this issue warrants further investigation. 

The authors wish to thank Martin Bucher for instruc- 
tive conversations. 
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